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ALGEBRAIC BOUNDARIES OF HILBERT'S SOS CONES 

GRIGORIY BLEKHERMAN, JONATHAN HAUENSTEIN, 

JOHN CHRISTIAN OTTEM, KRISTIAN RANESTAD 

AND BERND STURMFELS 

Abstract. We study the geometry underlying the difference between non- 
negative polynomials and sums of squares. The hypersurfaces that discrim- 
inate these two cones for ternary sextics and quaternary quartics are shown 
to be Noether-Lefschetz loci of K3 surfaces. The projective duals of these hy- 
persurfaces are defined by rank constraints on Hankel matrices. We compute 
their degrees using numerical algebraic geometry, thereby verifying results 
due to Maulik and Pandharipande. The non-SOS extreme rays of the two 



■^^ ■ cones of non-negative forms are parametrized respectively by the Severi va- 

(-H I riety of plane rational sextics and by the variety of quartic symmetroids 



1. Introduction 



^ ' A fundamental object in convex algebraic geometry is the cone S„ 2d of lio- 

^ . mogeneous polynomials of degree 2d in M.[xi, . . . , x„] that are sums of squares 

QQ ! (SOS). Hilbert [13j showed that the cones S3 g and S44 are strictly contained in 

the corresponding cones P^^q and P4^4 of non-negative polynomials. Blekherman 
t~^ ■ [5] furnished a geometric explanation for this containment. In spite of his recent 

progress, the geometry of the sets Pa^eV^a.e and P4^4\S4 4 remains mysterious. 
We here extend known results on Hilbert's SOS cones by characterizing their 
algebraic boundaries, that is, the hypersurfaces that arise as Zariski closures 
of their topological boundaries. The algebraic boundary of the cone Pn,2d of 
non- negative polynomials is the discriminant [I9], and this is also always one 
5^ I component in the algebraic boundary of S„^2d- The discriminant has degree 

n{2d — 1)"~^, which equals 75 for Ea^e and 108 for ^4^4. What we are interested 
in are the other components in the algebraic boundary of the SOS cones. 

Theorem 1. The algebraic boundary of S36 has a unique non-discriminant 
component. It has degree 83200 and consists of forms that are sums of three 
squares of cubics. Similarly, the algebraic boundary of ^4^4 has a unique non- 
discriminant component. It has degree 38475 and consists of forms that are 
sums of four squares of quadrics. Both hypersurfaces define Noether-Lefschetz 
divisors in moduli spaces of K3 surfaces. 
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Our characterization of these algebraic boundaries in terms of sums of few 
squares is a consequence of Corollaries 5.3 and 6.5]. What is new here is 
the connection to K3 surfaces, which elucidates the hypersurface of ternary 
sextics that are rank three quadrics in cubic forms, and the hypersurface of 
quartic forms in 4 variables that are rank four quadrics in quadratic forms. 
Their degrees are coefficients in the modular forms derived by Maulik and 
Pandharipande in their paper on Gromov-Witten and Noether-Lefschetz theory 
[IS] . In Section 2 we explain these concepts and present the proof of Theorem [TJ 

Section 3 is concerned with the cone dual to S„ 2d and with the dual varieties 
to our Noether-Lefschetz hypersurfaces in Theorem [H Each of them is a deter- 
minantal variety, defined by rank constraints on a 10 x 10-Hankel matrix, and it 
is parametrized by a Grassmannian via the global residue map in |6l §1.6]. We 
note that Hankel matrices are also known as moment matrices or catalecticants. 

Section 4 features another appearance of a Gromov-Witten number [iQ\ in 
convex algebraic geometry. Building on work of Reznick [23], we shall prove: 



Theorem 2. The Zariski closure of the set of extreme rays of P3fi\E3fi is the 
Severi variety of rational sextic curves in the projective plane P^. This Severi 
variety has dimension 17 and degree 26312976 in the P^"^ of all sextic curves. 



We also determine the analogous variety of extreme rays for quartics in P^: 

Theorem 3. The Zariski closure of the set of extreme rays of P4^4\S44 is 
the variety of quartic symmetroids in P^, that is, the surfaces whose defining 
polynomial is the determinant of a symmetric 4 x A-matrix of linear forms. This 
variety has dimension 24 in the P^^ of all quartic surfaces. 

Section 5 offers an experimental study of the objects in this paper using 
numerical algebraic geometry. We demonstrate that the degrees 83200 and 
38475 in Theorem [1] can be found from scratch using the software Bertini 
|4j. This provides computational validation for the cited results by Maulik and 
Pandharipande [18] . Motivated by Theorem [3l we also show how to compute a 
symmetric determinantal representation ([7]) for a given quartic symmetroid. 

A question one might ask is: What's the point of integers such as 38475? 
One answer is that the exact determination of such degrees signifies an under- 
standing of deep geometric structures that can be applied to a wider range of 
subsequent problems. A famous example is the number 3264 of plane conies 
that are tangent to five given conies. The finding of that particular integer 
in the 19th century led to the development of intersection theory in the 20th 
century, and ultimately to numerical algebraic geometry in the 21st century. 
To be more specific, our theorems above furnish novel geometric representa- 
tions of boundary sums of squares that are strictly positive, and of extremal 
non-negative polynomials that are not sums of squares. Apart from its intrinsic 
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appeal within algebraic geometry, we expect that our approach, with its focus 
on exphcit degrees, will be useful for applications in optimization and beyond. 

2. Noether-Lefschetz Loci of K3 Surfaces 

Every smooth quartic surface in P^ is a K3 surface. In our study of Hilbert's 
cone E4_4 we care about quartic surfaces containing an elliptic curve of degree 
4. As we shall see, these are the quartics that are sums of four squares. K3 
surfaces also arise as double covers of P^ ramified along a smooth sextic curve. 
In our study of Ss^g we care about K3 surfaces whose associated plane sextic 
is a sum of three squares. This constraint on K3 surfaces also appeared in the 
proof by Colliot-Thelene [8] that a general sextic in Es^g is a sum of four but 
not three squares of rational functions. 

Our point of departure is the Noether-Lefschetz Theorem (cf. [H]) which states 
that a general quartic surface 5* in P^ has Picard number 1. In particular, the 
classical result by Noether [20] and Lefschetz [T7] states that every irreducible 
curve on S is the intersection of S with another surface in P^. This has been 
extended to general polarized K3 surfaces, i.e. K3 surfaces 5* with an ampel 
divisor A. For each even I > 2, the moduli space Mi of K3 surfaces with a 
polarization A of degree A"^ = I has dimension 19 and is irreducible. For the 
general surface 5* in M;, the Picard group is generated by A. The locus in the 
moduli space Mi where the Picard number of S increases to 2 has codimension 
one. We are here interested in one irreducible component of that locus in 
M2, and also in M4. The relevant enumerative geometry was developed only 
recently, by Maulik and Pandharipande jTB], and our result rests on theirs. 



Proof of TheoremUi It was shown in [5] that SSaeV^Ps^e consists of ternary 
sextics F that are sums of three squares over R. Over the complex numbers C, 
such a sextic F is a rank three quadric in cubic forms, so it can be written as 

F = fh-g"^ where f,g, he C[xi,X2,X3\3. 

Let S be the surface of bidegree (2, 3) in P^xP^ defined by the polynomial 

G = fs^ + 2gst + ht^. 

If /, g and h are general, then the surface S is smooth. The canonical divisor on 
P^ X P^ has bidegree (—2, —3), so, by the adjunction formula, S* is a K3-surface. 
The projection S* — )■ P^ is two-to-one, ramified along the curve {F = 0} C P^. 
Up to the actions of SL{2, C) and SL{3, C), there is an 18-dimensional family 
of surfaces of bidegree (2,3) in P^ x P^. These surfaces determine a divisor 
D{2, 3) in the moduli space M2 of K3-surfaces with a polarization of degree 2. 
The Picard group of a general point S in this divisor has rank < 2. 

Let yl C S" be the preimage of a general line in P^, and let B he a. general fiber 
of the projection S* — ?■ P^. The classes of the curves A and B are independent in 
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the Picard group of S. Therefore the Picard group has rank at least 2 for every 
surface S G D{2,3). The curves A and B determine the intersection matrix 

( A^ A.B\ _ (2 3 

Conversely, any K3-surface with Picard group generated by classes A and 
B having the intersection matrix ([1]) has a natural embedding in P^ x P^ as a 
divisor of bidegree (2, 3): The linear system |A + i?| defines an embedding of S 
into the Segre variety P^ x P^ C P^ ([221 Proposition 7.15 and Example 7.19]). 

A general pencil of plane sextic curves contains a finite number of curves 
that are ramification loci of K3 double covers with Picard group of rank 2 and 
intersection matrix with a given discriminant. In our situation, this number is 
the degree of the hypersurface that forms the Zariski closure of S3 6\-P3,6- 

We shall derive this number from results of [18]. Let i? be a general surface 
of bidegree (2, 6) in P^ x P^, and let X be the double cover of P^ x P^ ramified 
along R. The general fiber of the projection X — )■ P^ is a curve 11 = P^ of 
K3 surfaces with a polarization A of degree I = A"^ = 2. Since the surface 
R has degree two in the first factor, the curve 11 defines a conic in the space 
of ternary sextics. Section 6 in [TS] computes the Noether-Lefschetz number 
A^L^3 of pairs {S, B) where [S] G 11, and B is the class of a curve of genus 
g{B) = 1 on S* and degree A-B = 3. By adjunction, the self-intersection equals 
B^ = 2g{B) — 2 = 0, and hence the intersection matrix is as above. 

There are two curve classes on the K3 surface S, namely B and 3A — B, 
that have self- intersection and intersection number A- B = A- {3 A — B) = 3. 
Therefore each such surface 5* appears twice in the count of ^3\- Also, the 
curve n is a conic in the space of plane sextic curves. So, to the get the count 
of surfaces S in the Noether-Lefschetz locus for a line in the space of sextics, 
we altogether must divide the number NL^^ by 4. 

In [ISl Corollary 3], the Noether-Lefschetz number NL^,^-.^^ is expressed 

as the coefficient of the monomial q^ in the expansion of a modular form 
0p of weight 21/2 as a power series in g^/^', where I = A^ is the degree of the 
polarization. The exponent of the relevant monomial is 6 = Ai{g{B),A-B)/2l, 

I A-B 



where A; ((7(5), A -5) is the discriminant of intersection matrix > ^ „ „ 

For the conic 11 in the space of sextics, the modular form 6^ has the expansion 

e" = -1 + 150g + 1248gt + lOSGOOg^ + 332800gt + 5113200^^ + . . . . 

In our case, we have 6 = 9/4, since / = 2 and the intersection matrix ([1]) has 
discriminant 9. We conclude that the number 

-A^L" = -332800 = 83200 

4 1.3 4 

equals the degree of the hypersurface of sextics that are sums of three squares. 
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We now come to the case of quartic surfaces in F^. It was shown in [5] that 
9S4 4\c}P4^4 consists of quartic forms F that are sums of four squares over M. 
Over the complex numbers C, such a quartic F is a rank 4 quadric in quadrics: 

(2) F = fg — hk = det H 1 ior some f,g,h,k e C[xi,X2,X3,X4\2. 

The K3 surface S defined by F contains two distinct pencils of elliptic curves on 
S, one defined by the rows and one by the columns of the 2x2 matrix. Up to 
the action of S'L(4,C), the determinantal quartics (|2]) form an 18-dimensional 
family, hence a divisor in the moduli space M4. A general surface S in this 
family has Picard rank 2, and its Picard group is generated by the class of a 
plane section and the class of an elliptic curve in one of the two elliptic pencils. 
Conversely, any smooth quartic surface 5* that contains an elliptic quartic 
curve is defined by a determinant F as in ([2]). This form of the equation is 
therefore characterized by the intersection matrix of S. Let A be the class of 
the plane section of S in P^ and let B and E be the classes of the curves in the 
two elliptic pencils. Then A and B = 2A — E have intersection numbers 

, . ( A^ A-B\ _ ( A^ A-E\ f A A 



\A-B B' J \A-E E^ J \A 

For general S, the classes A and B generate the Picard group and have inter- 
section matrix ([3]) with discriminant A4(l,4) = 16. Let 11 be a general linear 
pencil of quartic surfaces in P^. The Noether-Lefschetz number NL^^ counts 
pairs {S,B) where [S] G LI and i? is a curve class on S of degree 4 and genus 
g{B) = 1. Since there are two classes of such curves on S, we get the number 
of surfaces in the pencil containing such a curve class, by dividing NL^^ by 2. 
As above, the number NL^,^s ^^ is the coefficient of the monomial q^ in the 
expansion of a modular form Bp of weight 21/2 as a power series in g^/^', where 
I = A"^ is the degree of the polarization. Here 

^ A,{g{B),A-B) ^ 16 ^ 
8 8 

The modular form for the general line LI in the space of quartic surfaces equals 

ej = -l + 108g + 320gt +5016gt + 76950g2 + 136512g¥ + ... 



This was shown in [T8| Theorem 2]. We conclude that the degree of the hyper- 
surface of homogeneous quartics in 4 unknowns that are sums of 4 squares is 

-NL^, = -76950 = 38475. 
2 ^'^ 2 

This completes the proof of Theorem [H D 
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Remark 4. It was pointed out to us by Giorgio Ottaviani that the smooth 
ternary sextics that are rank three quadrics in cubic forms are known to co- 
incide with the smooth sextics that have an effective even theta characteristic 
(cf. [211 Proposition 8.4]). Thus the algebraic boundary of Hilbert's SOS cone 
for ternary sextics is also related to the theta locus in the moduli space A^iq. □ 

3. Rank Conditions on Hankel Matrices 

We now consider the convex cone (Sn,2(i)^ dual to the cone S„_2d- Its ele- 
ments are the linear forms £ on M.[xi, . . . , Xn]2d that are non-negative on squares. 
Each such linear form i is represented by its associated quadratic form on 
M[a;i, . . . , Xn\d, which is defined by / i— )■ ^(/^). The symmetric matrix which ex- 
presses this quadratic form with respect to the monomial basis of R[xi, . . . , Xn]d 
is denoted Hg, and it is called the Hankel matrix of i. It has format ("^^~ ) x 



' ) , and its rows and columns are indexed by elements of {{ii, i2, 



(4) He 



N^ : ii + i2 + ■ ■ ■ + in = d}. We shall examine the two cases of interest. 
The Hankel matrix for ternary sextics {n = d = 3) is the lOxlO-matrix 

O006 C^OIS '^024 ^033 C^IOS C^IU fll23 0-204 0,213 0303 

OOI5 O024 ^033 O042 Oll4 Ol23 O132 0213 0,222 ^312 

O024 O033 Oo42 O051 O123 (3132 Ol41 ^222 0231 O321 

O033 Oo42 O051 O06O ^132 O141 Oisq 0231 0240 O330 

Ol05 O114 O123 Oi32 O204 O213 O222 O303 O312 O402 

O114 O123 O132 Oi4i O213 O222 O231 O312 O321 O411 

O123 0,132 Oi4i Oi50 O222 O231 O240 O321 O330 O420 

O204 O213 O222 O231 O303 O312 O321 O402 O411 O501 

O213 O222 O23I O240 O312 O321 O330 O411 O420 O510 

O303 O312 O321 O330 O402 O411 O420 O501 O510 OgOO. 

The Hankel matrix for quaternary quartics {n = 4, d = 2) also has size 10x10: 

O0004 Oooi3 O0022 O0103 O0112 O0202 O1003 O1012 O1102 O2002 

O0OI3 O0022 O003I O0112 O012I O02II O1012 O102I Oiiii O2011 

O0022 O003I Oo040 O012I O0130 O022O O1021 O1030 O1120 O2020 

O0IO3 O0112 O012I O0202 O02II O0301 O1102 Oiiii O120I O2101 

O0II2 O012I O0130 O0211 O0220 O0310 Oiiii O1120 O1210 O2110 

O0202 O0211 O022O O030I O0310 Oo400 O1201 O1210 O1300 O22OO 

O1OO3 O1012 O1021 O1102 Oiiii O1201 O2002 O2011 O2101 O3001 

O1012 O102I O1030 Oiiii O1120 O1210 O20II O202O O2110 O3010 

O1IO2 Oiiii O1120 O120I O1210 O1300 O210I O2110 O220O O3100 

.O2OO2 O20II O202O O210I O2110 O220O O3001 O3010 O3100 04000. 

We note that what we call Hankel matrix is known as moment matrix in the lit- 
erature on optimization and functional analysis, and it is known as (symmetric) 
catalecticant in the literature on commutative algebra and algebraic geometry. 



(5) He 
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The dual cone (Ss^s)^ is the spectrahedron consisting of all positive semidefi- 
nite Hankel matrices dl]). The dual cone (£42)'^ is the spectrahedron consisting 
of all positive semidefinite matrices ([5]). This convex duality offers a way of 
representing our Noether-Lefschetz loci via their projective dual varieties. 

Theorem 5. The Hankel matrices ^ having rank < 7 constitute a rational 
projective variety of dimension 21 and degree 2640. Its dual is the hypersurface 
of sums of three squares of cubics. Likewise, the Hankel matrices ^ having 
rank < 6 constitute a rational projective variety of dimension 24 and degree 
28314. Its dual is the hypersurface of sums of four squares of quadrics. 

Proof. The fact that these varieties are rational and irreducible of the asserted 
dimensions can be seen as follows. Consider the Grassmannian Gr(3, 10) which 
parametrizes three-dimensional linear subspaces F of the 10-dimensional space 
M[a:i, X2, Xsjs of ternary cubics. This Grassmannian is rational and its dimension 
equals 21. The global residue in P^, as defined in [6l §1.6], specifies a rational 
map F n- Res(F) from Gr(3, 10) into F{(R[xi,X2,X3]q)*) ~ P^''. The base locus 
of this map is the resultant of three ternary cubics, so Res^i?) is well-defined 
whenever the ideal (F) is a complete intersection in ]R[xi,X2,X3]. The value 
Res(i?) (P) of this linear form on a ternary sextic P is the image of P modulo 
the ideal (F), and it can be computed via any Grobner basis normal form. Our 
map F I— )■ £ is birational because it has an explicit inverse: F = kernel (iif^). 
The inverse simply maps the rank 7 Hankel matrix representing i to its kernel. 
The situation is entirely analogous for n = A,d = 2. Here we consider the 24- 
dimensional Grassmannian Gr(4, 10) which parametrizes 4-dimensional linear 
subspaces F W[xi, X2, X3, X4]2. The global residue in P^ specifies a rational map 

Gr(4,10) — ^ P((M[xi,X2,X3,X4]4)*) ^P^^ F ^ Res(ir). 

This map is birational onto its image, the variety of rank 6 Hankel matrices 
([S]), and the inverse of that map takes a rank 6 Hankel matrix (^ to its kernel. 
To determine the degrees of our two Hankel determinantal varieties, we argue 
as follows. The variety Sr of all symmetric 10 x 10-matrices of rank < r is known 
to be irreducible and arithmetically Cohen-Macaulay, it has codimension ( ~^) , 
and its degree is given by the following formula due to Harris and Tu [10]: 

(a) --(-) ^^((/:i.)/r■r))■ 

Thus Sr has codimension 6 and degree 2640 for r = 7, and it has codimen- 
sion 10 and degree 28314 for r = 6. The projective linear subspace of Hankel 
matrices (jl]) has dimension 27. Its intersection with 5*7 was seen to have dimen- 
sion 21. Hence the intersection has the expected codimension 6 and is proper. 
That the intersection is proper ensures that the degree remains 2640. Likewise, 
the projective linear subspace of Hankel matrices ([5]) has dimension 34, and 
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its intersection with 5*6 has dimension 24. The intersection has the expected 
codimension 10, and we conclude as before that the degree equals 28314. 

It remains to be seen that the two Hankel determinantal varieties are projec- 
tively dual to the Noether-Lefschetz hypersurfaces in Theorem [1] This follows 
from [HI Corollary 5.2] for sextics curves in P^ and from [5i Corollary 5.7] for 
quartic surfaces in P^. These results characterize the relevant extreme rays of 
Eg g and S4 4 respectively. These rays are dual to the hyperplanes that support 
dTj^fi and c?E4^4 at smooth points representing strictly positive polynomials. By 
passing to the Zariski closures, we conclude that the algebraic boundaries of 
Eg^eV-Ps.e and E4 4\P4^4 are projectively dual to the Hankel determinantal vari- 
eties above. For a general introduction to the relationship between projective 
duality and cone duality in convex algebraic geometry we refer to [25]. D 

Remark 6. In the space P(Sym^V) of quadratic forms on a 10- dimensional 
vector space V*, the subvariety Sr of forms of rank < r is the dual variety to 
SlQ_r C P(Sym^y*). Identifying V with ternary cubics, the space of 10 x 10 
Hankel matrices 01]) form a 27-dimensional linear subspace H C P(Sym^y*). 
For r < 3, we have dim(5'r) < 27 and the variety dual to i/io-r = SIq_^. fl H 
coincides with the image S,. of the birational projection of 5*^ into H*. That 
image is the variety of sextics that are quadrics of rank < r in cubics. Further- 
more, when r < 2 the projection from S* to S^ is a morphism, so the degrees of 
these two varieties coincide. When r = 3, the projection is not a morphism and 
the degree drops to 83200. A similar analysis works for V = ]R[xi,X2,X3,X4]2 
with r < 4. D 

4. Extreme Non-negative Forms 

For each of Hilbert's two critical cases, in Section 2 we examined the hyper- 
surface separating T,n,2d and Pn,2d\^n,2d- In this section we take an alternative 
look at this separation, namely, we focus on the extreme rays of the cone Pn,2d 
of non-negative forms that do not lie in the SOS subcone S„ 2d- We begin with 
the following result on zeros of non-negative forms in the two Hilbert cases. 

Proposition 7. Let p be a non-negative form in P^q or P^^^. If p has more 
than 10 zeros, then p has infinitely many zeros and it is a sum of squares. 

Proof. The statement for Pgg was proved by Choi, Lam and Reznick in [7|. 
They also showed the statement for the cone P4^4 but with "11 zeros" instead 
of "10 zeros" . To reduce the number from 11 to 10, we use Kharlamov's theorem 
in [15] which states that the number of connected components of any quartic 
surface in real projective 3-space is < 10. See also Rohn's classical work [23]. D 

Recall that a face of a closed convex set i^ in a finite-dimensional real vector 
space is exposed if it is the intersection of K with a supporting hyperplane. The 
extreme rays of K lie in the closure (and hence in the Zariski closure) of the 
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set of exposed extreme rays [27] . A polynomial p G Pn,2d\^n,2d that generates 
an extreme exposed ray of Pn,2d will be called an extreme non-negative form. 

Our first goal is to prove Theorem [2], which characterizes the Zariski closure 
of the semi- algebraic set of all extreme non- negative forms for n = d = 3. 

Proof of Theorem\^ Suppose p G P3,6\S3,6 is an extreme form. By [231, Lemma 
7.1], the polynomial p is irreducible. Moreover, we claim that |Vr(p)| > 10. It 
is not hard to show that p is an extreme non- negative form if and only if Vm(p) 
is maximal among all forms in Pn,2d- In other words, if p is an extreme non- 
negative form and Vr(p) C V^{q) for some q G Pn,2d then q = Xp for some 
A G M. Now suppose that |Vr(p)| < 9. Then there is a ternary cubic q G Pa^s 
that vanishes on Vr(p). We have q^ G Ps^e and Vr(p) C Vu^q). This contradicts 
maximality of Vr{p). By Proposition [7] we conclude that |Vr(p)| = 10. 

Let C be the sextic curve in the complex projective plane P^ defined by p = 0. 
Since C is irreducible, it must have non- negative genus. Each point in Vr(p) 
is a singular point of the complex curve C. As this gives C ten singularities, 
it follows by the genus formula that C can have no more singularities, and 
furthermore that all of the real zeros of p are ordinary singularities. The genus 
of C is zero and therefore it is an irreducible rational curve. 

Let Sgfi denote the Severi variety of rational sextic curves in P^. We have 
shown that iS6,o contains the semi-algebraic set of extreme forms in Ps^gV^s^e- 
This is a subvariety in the P^^ of ternary sextics. The Severi variety iS6,o is 
known to be irreducible, and the general member C has exactly 10 nodes. 
Moreover, that set of 10 nodes in P^ uniquely identifies the rational curve C. 

Each rational sextic curve in P^ is the image of a morphism P^ — )■ P^ defined 
by three binary forms of degree 6. To choose these, we have 3-7 = 21 degrees 
of freedom. However, the image in P^ is invariant under the natural action 
of the 4-dimensional group GL{2, C) on the parametrization, and hence Sq^ 
has dimension 21 — 4 = 17. The degree of 56,0 is the number of rational 
sextics passing through 17 given points in P^, which is one of the Gromov- 
Witten numbers of P^. For rational curves, these numbers were computed by 
Kontsevich and Manin |T6] using an explicit recursion formula equivalent to the 
WDVV equations. From their recursion, one gets degree(iS6,o) = 26312976. 

To complete the proof, it remains to be shown that the semi-algebraic set 
of extreme forms in P3,6\S3,6 is Zariski dense in the Severi variety Sq^. We 
deduce this from j23l Theorem 4.1 and Section 5]. There, starting with a 
specific set F of 8 points in P^, an explicit 1-parameter family of non- negative 
sextics with 10 zeros, 8 of which come from F, was constructed using Hilbert's 
Method. Furthermore, by Theorem 4.1, Hilbert's Method can be applied to any 
8 point configuration in the neighborhood of F. By a continuity argument, all 
8 point configurations sufficiently close to F will also have a 1-parameter family 
of non-negative forms with 10 zeros. All such forms are exposed extreme rays. 
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This identifies a semi-algebraic set of extreme non-negative forms having 
dimension 16 + 1 = 17. We conclude that this set is Zariski dense in Sq^. D 

Remark 8. Our analysis implies the following result concerning SPs^gV^s.e- 
All exposed extreme rays are sextics with ten acnodes, and all extreme rays 
are limits of sextics with ten acnodes. This proves the second part of Reznick's 
Conjecture 7.9 in [23]. Indeed, in the second paragraph of the above proof we 
saw that C has ten ordinary singularities. These cannot be cusps since p > 0. 
Hence they have to be what is classically called acnodes, or round zeros in 



Our next goal is to derive Theorem [3l the analogue to Theorem [2] for quartic 
surfaces in P^. The role of the Severi variety ^6,0 is now played by the variety 
QS of quartic symmetroids, i. e. the surfaces whose defining polynomial equals 

(7) F{xi,X2,X3,X4) = det{AiXi + A2X2 + A3X3 + A4X4), 

where Ai, A2, A3, A^ are complex symmetric 4 x 4- matrices. 

Lemma 9. The variety QS is irreducible and has codimension 10 in P^^. 

Proof. Each of the four symmetric matrices Ai has 10 free parameters. The 
formula ([7]) expresses the 35 coefficients of F as quartic polynomials in the 40 
parameters, and hence defines a rational map P^^ --■> P^^. Our variety QS is 
the Zariski closure of the image of this map, and so it is irreducible. To compute 
its dimension, we form the 35 x 40 Jacobian matrix of the parametrization. By 
evaluating at a generic point {Ai, . . . , A4), we find that the Jacobian matrix 
has rank 25. Hence the dimension of the symmetroid variety QS C P^^ is 24. 
For a theoretical argument see [HI page 168, Chapter IX. 101]. D 

A general complex symmetroid S has 10 nodes, but not every 10- nodal quartic 
in P^ is a symmetroid. To identify symmetroids, we employ the following lemma 
from Jessop's classical treatise [2] on singular quartic surfaces. Let S" be a 10- 
nodal quartic with a node at p = (0 : : : 1). Its defining polynomial equals 
F = fx\ + 2gx4 + h where f,g,h G C[xi,X2,a;3] are homogeneous of degrees 
2,3,4 respectively. The projection of 5* from p is a double cover of the plane 
with coordinates xi, X2, X3 ramified along the sextic curve Cp defined by g^ — fh. 
The curve Cp has nodes exactly at the image of the nodes on S that are distinct 
from p. Since no three nodes on S are collinear, the curve Cp has 9 nodes in 
P^. The following result appears on page 14 in Chapter 1.8 of 



Lemma 10. // the sextic ramification curve Cp is the union of two smooth 
cubics that intersect in 9 distinct points, then the quartic surface S is a sym- 
metroid and, moreover, the ramification curve Cq for the projection from any 
node q on S is the union of two smooth cubic curves. 
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Proof of Theorem\^ Let £ denote the seinialgebraic set of all non-negative ex- 
treme forms F in P4^4\S4 4. Each F E S satisfies |Vm(F)| = 10, by Proposition 
[7] and the same argument as in the first paragraph in the proof of Theorem [2j 
Thus S consists of those real quartic surfaces in P^ that have precisely 10 real 
points. 

We shall prove that £^ is a subset of QS. Let F E £ and S = Vc{F) the 
corresponding complex surface. Then S" is a real quartic with 10 nodes, and 
these nodes are real. Our goal is to show that S" is a symmetroid over C If 
p G Pg is one of the nodes of F, then the ramification curve Cp is a real sextic 
curve with 9 real nodes at the image of the nodes distinct from p. Since the 
nodes on S are the only real points, these nodes are the only real points on Cp. 

Through any nine of the nodes of S there is a real quadratic surface. This 
quadric is unique; otherwise there is a real quadric through all ten nodes and 
F is not extreme. Let g be a node on S distinct from p and A a real quadratic 
form vanishing on all nodes on 5* except q. Consider the pencil of quartic forms 

Ft = F + tA^ for t G M. 

Suppose p = (0 : : : 1). The polynomial A has the form ux^ + v, where 
u,v E M[a;i, a;2, x^] have degree 1 and 2. The equation of F^ is then given by 

Ft = {f + tu'^)xl + 2{g + tuv)x4 + h + tv^. 

Any surface St = {Ft = 0} has at least 9 real singular points, namely the nodes 
of S other than q. Since F is non-negative. Ft is non-negative for t > with 
zeros precisely at the 9 nodes. On the other hand, F has an additional zero at 
q. Since A"^ is positive at q, the real surface {Ft = 0} must have a 2-dimensional 
component when t < 0. Projecting from p we get a pencil of ramification loci 
Cp(t). In the above notation, this family of sextic curves is defined by the forms 

Gt = fh — g'^ + t{hu'^ — 2guv + fv'^) E M.[xi,X2,X3]q. 

The curves in this pencil have common nodes at eight real points pi, . . . ,ps in 
the plane P^, namely the images of the nodes on S other than p and q. 

Consider the vector space V of real sextic forms that are singular at pi, . . . ,p8- 
Since each pi imposes 3 linear conditions, we have dim V^ > 28 — 3 ■ 8 = 4. We 
claim that dim V = 4. To see this, consider a general curve Cp{t) with t > 0. 
It has only eight real points, so as a complex curve it is irreducible and smooth 
outside the eight nodes. Hence the geometric genus of Cp(t) is 2. Let X denote 
the blow-up of the plane in the points pi, . . . ,ps, and denote by C the strict 
transform of Cp{t) on X. By Riemann-Roch, diia H^(Ox{C)\„) = 3, since 
C^ = 4. Combined with the cohomology of the exact sequence 

we conclude that dim V = dimif°(X, Ox{C)) < 4, and hence dimy = 4. 
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The pencil Mj/ci, k2} of real cubic forms through the eight points pi, . . . ,ps 
determine a 3-dimensional subspace U = Mjfc^, kik2, k"^} of V, while the sextic 
forms Gt span a 2-dimensional subspace L of V. Since Gt has no real zeros 
except the nodes when t > 0, we see that L is not contained in U. Hence L 
and U intersect in a 1-dimensional subspace of V, so there exists a unique value 
io e IR such that Cp(to) = Ki ■ K2, where Ki,K2 G C[xi,X2,X3]3. 

We now have two possibilities: either Ki and K2 are both real, or Ki and 
K2 are complex conjugates. We claim that the latter is the case. Consider the 
intersection {Gt = 0} r\ {Ki ■ K2 = 0}. This scheme is the union of a scheme of 
length 32 supported on the 8 nodes and a scheme Z of length 6 ■ 6 — 4 ■ 8 = 4. Its 
defining ideal {fh — g"^, hu^ — 2guv + fv"^) contains the square {gu — 2/f )^, and 
thus each of its points has even length. Hence Z is either one point of length 
4 or two points of length 2. Since the general Gt does not contain the ninth 
intersection point of Ki and i^2, each component of Z is contained in only one 
of the Ki. In particular, since Ki fl Cp{t) contains a scheme of length 2 disjoint 
from the points pi, . . . ,P8) this shows that Z has two points, one in each of the 
Ki. If both Zi were real then KiH Z would be real, contradicting the fact that 
Gt has only 8 real points. We conclude that the two cubics Ki^ K2 are complex 
conjugates and their only real points are the 9 common intersection points. 

We now claim that to = 0. Indeed, if to < then St^ has 2-dimensional 
real components and the real points in the ramification locus Cp(to) would have 
dimension 1. If to > then St^ has only 9 real points and Cp{to) has only 8 real 
points. Since Cp(to) = Ki ■ K2 has 9 real zeros, it follows that to = 0. Using 
Jessop's Lemma [TOl we now conclude that F = Fq is a symmetroid. 

We have shown that the semi-algebraic set S is contained in the symmetroid 
variety QS. It remains to be proved that S is Zariski dense in QS. To see this, 
we start with any particular extreme quartic. For instance, take the following 
extreme quartic due to Choi, Lam and Reznick [71 Proof of Proposition 4.13]: 

(8) Fb = ^2^^^^j + b'^xfxjXk + {Ab'^-Ab-2)xiX2X3X4 for 1 < 6 < 2, 

where the sums are taken over all distinct pairs and triples of indices. The 
complex surface defined by Fi, has 10 nodes, namely, the points in Vm.{Fj,). Our 
proof above shows that F^ is a symmetroid. Since the Hessian of Ff, is positive 
definite at each of the 10 real points, we can now perturb these freely in a small 
open neighborhood inside the variety of 10-tuples of real points that are nodes of 
a symmetroid. Each corresponding quartic is real, non-negative and extreme. 
This adaptation of "Hilbert's method" constructs a semi- algebraic family of 
dimension 24 in £^. We conclude that QS is the Zariski closure oi S. D 

Our proof raises the question whether Lemma [TOl can be turned into an algo- 
rithm. To be precise, given an extreme quartic, such as ([8]), what is a practical 
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method for computing a complex symmetric determinantal representation ([7])? 
We shall address this question in the second half of the next section. 

5. Numerical Algebraic Geometry 

We verified the results of Theorems [Hand [5] using the algorithms implemented 
in Bertini |3]. In what follows we shall explain our methodology and findings. 
An introduction to numerical algebraic geometry can be found in |29] . 

The main computational method used in Bertini is homotopy continuation. 
Given a polynomial system F with the same number of variables and equations, 
basic homotopy continuation computes a finite set S of complex roots of F 
which contains the set of isolated roots. By "computes iS" we mean a numerical 
approximation of each point in S together with an algorithm for computing each 
point in S to arbitrary accuracy. The basic idea is to consider a parameterized 
family T of polynomial systems which contains F. One first computes the 
isolated roots of a sufficiently general member of J-", say G, and then tracks the 
solution paths starting with the isolated roots of G at t = 1 of the homotopy 

H{x,t) = F{x){l-t)+tG{x). 

The solution paths are tracked numerically using predictor-corrector methods. 
For enhanced numerical reliability, the adaptive step size and adaptive precision 
path tracking methods of [3] is used. The endpoints at t = of these paths 
can be computed to arbitrary accuracy using endgames with the set of finite 
endpoints being the set S. If F has finitely many roots, then S is the set of all 
roots of F. If the variety of F is not zero-dimensional, then the set of isolated 
roots of F is obtained from S using the local dimension test of [2]. 

Our computations to numerically verify the degrees in Theorem [1] only used 
basic homotopy continuation. For the Ss^g case, we computed the intersection 
of the set of rank three quadrics in cubics with a random line in the space P^^ 
of ternary sextics. In particular, for random p, g G C[xo, xi, a;2]6, we computed 
the complex values of s such that there exists f,g,hE C[a;o, xi, X2J3 with 

fh-g^ = p + sq. 

We used the two degrees of freedom in the parametrization of a rank three 
quadric in cubics by taking the coefficient of Xg in g and XqXi in / to be zero, and 
we dehomogenized by taking the coefficient of Xq in / to be one. The resulting 
system F = consists of 26 quadratic and two linear equations in 28 variables. 
Since the solution set of F = is invariant under the action of negating g, 
we considered F as a member of the family J-" of all polynomial systems in 
28 variables consisting of two linear and 26 quadratic polynomials which are 
invariant under this action. It is easy to verify that a general member of J-' has 
226 j-Qots, which consist of 2^^ orbits of order 2 under the action of negating 
g. We took the system G to be a dense linear product polynomial system [30] 
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with random coefficients which respected this action. By tracking one path 
from each of the 2^^ orbits, which took about 40 hours using 80 processors, this 
yielded 166400 points which correspond to 83200 distinct values of s. 

The S4 4 case of Theorem [1] was solved similarly, and the number 38475 was 
verified. We took advantage of the bi-homogeneous structure of the system 

fg — hk = p + sq. 

Numerical algebraic geometry can be used to compute all irreducible compo- 
nents of a complex algebraic variety. Here the methods combine the ability to 
compute isolated solutions with the use of random hyperplane sections. Each 
irreducible component ^ of F = is represented by a witness set which is a 
triple [F, C, W) where £ is a system of dim V random linear polynomials and 
W is the finite set consisting of the points of intersection of V with £ = 0. 

Briefly, the basic approach to compute a witness set for the irreducible com- 
ponents of -F = of dimension k is to ffist compute the isolated solutions W of 
F = Ck = where C^ is a system of k random linear polynomials. The set W 
is then partitioned into sets, each of which corresponds to the intersection of 
£fc = with an irreducible component of F = of dimension k. The cascade 
[28] and regenerative cascade [12] algorithms use a sequence of homotopies to 
compute the isolated solutions oi F = Ck = for all relevant values of k. 

We applied these techniques to verify the results of Theorem |5] concerning our 
10 X 10 Hankel matrices. Our computations combined the regenerative cascade 
algorithm with the numerical rank- deficiency method oi [I]. In short, if A{x) is 
an n X A^ matrix with polynomial entries, consider the polynomial system 



Fr = A{x) ■ B 



iN-r 



where B G C^^^ is random, I^-r is the {N — r) x (N — r) identity matrix, 
and S is an r X (A^ — r) matrix of unknowns. One computes the irreducible 
components of F^ = whose general fiber under the projection [x, S) 1— )■ x is 
zero-dimensional. The images of these components are the components of 

Sr{A) = {x : TankA{x)<r}. 

The degree of such degeneracy loci is then computed using the method of [TTj. 
The results on degree and codimension in Theorem [S] were thus verified, with 
the workhorse being the regenerative cascade algorithm. For instance, we ran 
Bertini for 12 hours on 80 processors to find that the variety of Hankel matrices 
(jlj) of rank < 7 is indeed irreducible of dimension 21 and degree 2640. 

We now shift gears and discuss the problem that arose at the end of Section 
4, namely, how to compute a symmetric determinantal representation ([7]) for 
a given extremal quartic F & £ d 9P4,4\E4^4. For a concrete example let 
us consider the Choi-Lam-Reznick quartic in ([8]) with h = 3/2. We found 
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that F3/2 = det(M)/7, where 7 = -54874315598400(735^ + 2201), with uj 
|\/— 10, and M is the symmetric matrix with entries 



mil = (-11844a; + 8100)xi + (3024a; + 13140)a;3, 

mi2 = (7980a; + 14820)a;3, 

mis = (19971a; - 17460)xi + (4494a; + 9600)a;3, 

mi4 = (-1596a; - 26790)x3 + (15561a; - 6840)x4, 

m22 = (30324a; - 7220)a;2 + (20216a; + 21660)^3, 

^' m23 = (20216a; + 21660)x2 + (6384a; + 27740)x3, 

m24 = (-20216a; - 21660)a;2 - 39710x3 + (7581a; - 21660)a;4, 

m33 = (-13230a; + 31860)xi + 39710x2 + (-28910a; + 29910)x3, 

m34 = -39710x2 + (25004a; -17100)x3 + ((5187/2)a;-1140)x4, 

m44 = 39710x2 + (-20216a; + 37905)x3 + (-30324a; + 27075)x4. 

A naive approach to obtaining such representations is to extend the nu- 
merical techniques introduced for quartic curves in ^F2\ §2]: after changing 
coordinates so that x\ appears with coefficient 1 in F, one assumes that Ai is 
the identity matrix, A2 an unknown diagonal matrix, and ^3 and ^4 arbitrary 
symmetric 4x4- matrices with unknown entries. The total number of unknowns 
is 4+10+10 = 24, so it matches the dimension of the symmetroid variety QS. 
With this, the identity ([7]) translates into a system of 34 polynomial equations 
in 24 unknowns. Solving these equations directly using Bertini is currently not 
possible. Since the system is overdetermined, Bertini actually uses a random 
subsystem which has a total degree of 3^4^^. The randomization destroys much 
of the underlying structure and solving this system is currently infeasible. 

In what follows, we outline a better algorithm based on the underlying ge- 
ometry of the problem. The input is a 10- nodal quartic surface S = {F = 0}. 
After changing coordinates, so that p = (0:0:0:1) is one of the nodes, the quartic 
has the form: 

F = fx1 + 2gx4, + h where /, (7, /i G M[xi,X2,X3]. 

The projection from p defines a double cover tt : S* — )■ P^ and the ramification 
locus is the sextic curve whose defining polynomial is fh — g^ and splits into 
a product of two complex conjugate cubic forms Ki,K2. The intersection of 
S with {Ki = 0}, regarded as a cubic cone in P^, is supported on the branch 
locus of the double cover and therefore equals two times a curve C of degree 6. 
The curve C has a triple point at the vertex p, its arithmetic genus is 3, and it 
is arithmetically Cohen- Macaulay. By the Hilbert-Burch Theorem, the ideal of 
C is generated by the 3 x 3- minors ^fi, . . . , 5(4 of a 3 x 4 matrix whose entries 
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are linear forms in C[xi,X2,X3,a;4]i: 
(10) 



^11 lu ha hi 

hi ^22 ^23 hi 
hi hi ^33 ^34 



The rows of this matrix give three hnear syzygies between the four cubics gi. 
Furthermore, F itself is in the ideal generated by these cubics, and so there 
is a hnear relation F = hgi + ■ ■ ■ + hg4- Hence the quartic F is equal, up to 
multiplication by a non-zero scalar in C, to the determinant of the matrix 

^1 —h h ~h 

hi hi hs hi 

hi hi hs hi 

hi hi hs hi 

To find a symmetric matrix M with the same property, we solve the linear 
system PL = (PL)^ for some matrix P G GL(4, C) and define M = PL. 

A numerical version of the above algorithm is almost exactly as explained 
above except that a basis for the ideal Ic of the genus 3 curve C is found by 
computing a large sample of points in the intersection {F = Ki = 0}, and then 
computing a basis gi, . . . ,g4 for the 4-dimensional space of cubic forms vanishing 
on this set. Next, a basis for the 3-dimensional set of linear syzygies between 
these cubics is computed. This yields the matrix in flTOj) whose 3x3 minors 
are the four cubics gi. For the quartic ([8]) with b = 3/2, we used Bertini 
to compute 100 random points in this intersection and then used standard 
numerical linear algebra algorithms. In total, it took 30 seconds to compute a 
symmetric determinantal representation for -F3/2. To four digits, with i = \/—l, 
the output we found is the symmetric matrix M with entries 



11 — 


(15.5378 + 5.6547i)i!i 


- (20.4008 - 


- 5.8116i)a;2 - 


(23.1956 + 16.9236i)a;3 + (12.4987 + 26.8206i)2!4, 


12 — 


(18.3458 - 5.8125i)a:i 


- (14.0867 - 


- 25.1505i)a:2 


- (35.0029 - 5.29481)2:3 + (36.1417 + 15.7167i)a:4 , 


13 — 


(11.6232 + 5.6624i)2^i 


- (15.6076 - 


- 5.9393i)a:2 - 


(17.3057 + 12.3685i)a;3 + (11.0079 + 22.8305i)i!4, 


14 — 


(25.7222 + 1.2098i)a;i - 


- (27.4233 - 


22.3864i)a:2 - 


- (45.8046 + 14.1068i)a;3 + (35.3836 + 37.8454i)a;4 


22 = 


(12.6315 - 18.4638i)a;i 


+ (9.6932 + 37.6953i)a;2 - 


- (26.0269 - 34.9909i)2:3 + (49.9098 - 16.2993i)a;4 


23 = 


(14.6285 - 3.0705i)a;i 


- (9.5983 - 


20.6203i)i2 - 


- (25.8489 - 4.2265i)2^3 + (31.1616 + 13.0794i)2:4, 


24 — 


(24.1544 - 17.3589i)a;i 


- (5.2755 - 


- 47.6528i)a:2 


- (52.3363 - 27.5281i)a:3 + (68.7313 + 6.4353i)2;4, 


33 = 


{ 8.5030 + 5.3275i)a:i 


- (11.9127 


- 5.6822i)a:2 


- (12.9473 + 8.9555i)a:3 + (9.6646 + 19.4288i)a:4, 


34 — 


(19.6130 + 2.9165i)i^l 


- (20.0754 - 


- 19.3371i)a^2 


- (34.2042 + 9.99111)2:3 + (30.7454 + 32.0581i)i:4, 


44 — 


(37.6831 - 10.7034i)a;i 


- (27.3051 


- 52.2852i)i^2 


- (80.4558 - 2.6947i)a;3 + (79.5452 + 43.7001i)i:4 



(11) 



The symbolic solution (^ and the numerical solution fITT]) are in the same 
equivalence class of symmetric matrix representations. In fact, we close with 
the result that the output of the algorithm is essentially unique, independant 
of the choice of node p and cubic form Kf. 

Proposition 11. For any 10-nodal symmetroid F G QS, the representation ^ 
is unique up to the natural action of GL(4, C) via Ai \-^ UAiU'^ fori = 1, 2, 3, 4. 
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Proof. Let M = Y^XiAi be a symmetric matrix such that F = det(M). Any 
three of the four rows of M determine a curve C by taking 3x3 minors. This 
gives a 4-dimensional hnear system Lm of curves of arithmetic genus 3 and 
degree 6 on S. The doubhng of any curve in Lm is the complete intersection of 
S and a cubic surface defined by a 3 x 3-symmetric submatrix of M. Conversely, 
the linear system determines the matrix M up to a change of basis. 

Each curve in Lm passes through all the nodes of S, and these are the common 
zeros of the curves in Lm- If S is the smooth K3 surface obtained by resolving 
the nodes, then by Riemann-Roch, Lm defines a complete linear system on S. 
Since Pic(S') is torsion-free, we see that Lm is uniquely determined as the linear 
system of degree 6 curves on S passing through all nodes and whose doubling 
form a complete intersection. Therefore the equivalence class of the symmetric 
matrix representation is also unique. D 
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